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The diagrammatic cumulant expansion for the periodic Anderson model with 
infinite Coulomb repulsion {U = oo) is considered here for an hypercubic lattice 
of infinite dimension (d — cxd). The nearest neighbor hopping of the uncorrelated 
electrons is described exactly by a conduction band, while two different models 
of hybridization are treated as a perturbation. The same type of simplifications 
obtained by Metzner for the cumulant expansion of the Hubbard model in the 
limit of d = oo, are shown to be also valid for the periodic Anderson model. The 
derivation of these properties had to be modified because of the exact treatment of 
the conduction band. 
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I. INTRODUCTION 

The periodic Anderson model (PAM) gives a schematic description of very important sj^stems 
with strongly correlated electrons, and there are several recent reviews on this subjectEJ. The 
model consists of a lattice with two localized electronic states at each site, strongly correlated by a 
Coulomb repulsion U, plus a band of uncorrelated conduction electrons (c-electrons) that hybridize 
with the localized electrons (f-electrons). The cumulant expansion of the PAM discussed in this 
work is a perturbative expansion around the atomic limit, and is an extensionB of the method 
originally employed by HubbardS to study his well known model of correlated itinerant electronsQ. 
He introduced operators Xjj,a =\ j,b){j,a |, that transform the local state | a) at site j into the 
local state | b) at the same site, and developed a diagrammatic method that circumvents the fact 
that the X operators are neither fermions nor bosons. |-|His work was the first application of the 
cumulant expansion to a quantum system of fermionsO and has several desirable properties: it 
seems to be a natural extension of the usual diagrammatic techniques, it does not have excluded 
site problems in the lattice summations and it is possible to derive a linked cluster expansion for 
the grand canonical potential( see refjj for more details). At the same time, the use of X operators 
allows to work in a subspace of the whole space of local states, that contains only the states that 
are relevant to the problem of interest. It seems therefore useful to understand better this method, 
and the main purpose of the present paper is to study the properties of the cumulant expansion 
for the PAM in the limit ai infinite spatial dimensions {d oo). 

Metzner and VoUhardttl have recently applied the limit d — s- cx) to strongly correlated fermion 
systems, providing non-trivial models of the Hubbard type that are substantially simpler to analyze. 
This development has stimulated many new works, and we could mentioned the very successful 
Dynamical Mean Field TheoryQ, that uses the local properties of the self energy for_jcL=r> oo as a 
starting point. Several methods have been employed to study the PAM_in this limitETEJ, and the 
cumulant expansion of the Hubbard model was reconsidered by Metzneiii3, showing that there is an 
important reduction in the number and type of cumulant diagrams that appear in the expansion 
when d — > oo in an hypercubic lattice. In the Hubbard model there is no hybridization and 
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the hopping Hamiltonian is used as perturbation, while the PAM is usually employed to model 
systems with a bandwidth much larger than the hybridization energies, and it is then preferable 
to diagonalize the conduction band and use the hybridization Hamiltonian as perturbation. It was 
then necessary to make substantial modifications to Metzner's derivation in order to apply it to 
the PAM when oo. 

In the following sections two types of hybridization models are considered: the purely "local" 
hybridization and the "nearest neighbor" (n.n.) one. In the n.n. hybridization there are electronic 
transitions between localized ("f") states at a site and conduction Wannier states at the nearest 
neighbor sites, while in the local hybridization those transitions only occur at the same site. The 
same type of diagram cancellation obtained by Metzner for the Hubbard model is obtained for 
these two hybridization models in the PAM. 

As the local hybridization between / and d or s electrons should vanish when the site has 
inversion symmetry, the n.n. hybridization discussed in this paper is more realistic than the local 
one, because it does not necessarily vanish for that type of site. The local approximation is usually 
employed because it makes the calculation simpler; and the changes introduced by pi interatomic 
hybridization have been already discussed in the case of a diatomic-molecule modeluJ. 



II. THE HYPERCUBIC LATTICE 



Let us coasider free electrons in an hypercubic lattice of dimension d, described by a 
HamiltonianU 

He = ^ ^ tn iYi C^^^ Cm,cr (2-1) 
n,m,(7 

with translational invariance. The position of the site n = (ni,...,nd) is given by R„ = 
X]q=i ^a, where \ea\ = a is the lattice parameter and the rij are integers. Only n.n. hop- 
ping shall be considered, so that in,m is non-zero and equal to —t when the components of m 
satisfy TOq = ria + Sa,-y for 7 = 1, 2, . . . , d. Employing the Wannier transformation 



C„,. = Ck,. exp[ik.R„] , (2.2) 



N 1. 

with k satisfying cyclic boundary conditions, one obtains 

Hc = Y.EuCl,Ck^. , (2.3) 

where 

d 

Ek = -2t^cos{ks a) . (2.4) 

s=l 

The electronic Green's functions (GF) for imaginary time are defined bjJ 

G.(k,r)^/(Ck,.(r) \ , (2.5) 



where Ck,o-(T) = exp[HcT] Ck.o- exp[—HcT] and the subindex -I- in Eq.(^.5|) is the usual chronolog- 
ical ordering of Fermi operators, with r increasing to the left. To emphasize that the hybridization 
was not considered in these GF, they will be deppted with G°(k,T), and after Fourier transforma- 
tion in the imaginary time r they are given byO G^{k,uj^) = —{itOu — i?k)~^- To obtain a finite 
average energy per site when d — > 00, it is necessary to renormalize -ffc, taking a finite t in the 
non-zero ii™ = —t/^/2d (i.e. when n and m are n.n.). In that limit the energy density of states 
is given byll^ 



Mi^) = ^exp 



p2 

2¥ ' 
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and 



{E{-k)f)^t^ . (2.7) 

The exact solution of He is used in the present treatment of the PAM, employing the hybridiza- 
tion as perturbation. To study the behavior of the cumulant expansion in the limit d oo, it is 
convenient to consider the direct space GF, given by 

G°(R„,„„w,)=^^Gl}(k,c^,)exp[ik.R„,^] , (2.8) 

k 

where R n.m — I^n R'mi as well as its dependence with the minimum number p of n.n. jumps 
necessary to go from m to n. When the GF are obtained for large d using th e ho pping as a 
perturbation, one immediately shows that (Rn.m? ^i/) — O|0^|, where 9 = t/V^d . A similar 



result is obtained for the "exact" G°(R„_m,^i/) of Eq. (2.8), but the discussion will be postponed 



to a latter section, in which a more general property will be shown. 



III. THE HAMILTONIAN OF THE PAM 



The Hamiltonian of the whole system is 

H = H, + Hf + Hh , (3.1) 

where He is the Hamiltonian of the conduction electrons (c-electrons) discussed in Section ||. The 
second term Hf = "^Ij a E^Xj^aa describes independent localized electrons (f-electrons), where 
a simple index j has been used to indicate the sites n. The last term is the hybridization 
Hamiltoniaru, giving the interaction between the c-electrons and the f-electrons: 

Hh=Y. (^J-^k,. C-k,. + V;^^, Cl^ X,,oa) . (3.2) 

The state space of the f-electrons, at a given site j, is spanned by four states: the vacuum state 
I j, 0), the two states] j, cr) of one f-electron with spin component a and the state | j, 2) with two 
electrons of opposite spin. In the limit of infinite electronic repulsion {U — oo) considered here, the 
state I j, 2) is always empty, and can be projected out of the space. The operator Xj^Q^ destroys 
the electron in state | j, a) leaving the site in the vacuum state | j, 0) with energy i?j,o — 0; ^^nd its 
Hermitian conjugate ^joo- — ^i,<TO reverses that process. The Hubbard operators Xj^ab are not 
usually Fermi or Bose operators, and a product rule: 

A'j.ab Xj.cd — Sbc Xj,ad , (3-3) 

should be employed when the two operators are at the same site. For different sites one should first 
classify the operators in two families: Xj^ab is of the "Fermi type" when | a) and | b) differ by an odd 
number of Fermions, and it is of the "Bose type" when they differ by an even number of Fermions. 
At different sites, two X-operators of the Fermi type anticommute, and they commute when at 
least one of them is of the Bose type. To complete the definition, the Xj^ab should anticommute 
(commute) with the G^ ^ and Gk,o- when they are of the Fermi type (Bose type). Sometimes it 
will be convenient to employ ^(7) to describe the Hubbard operators, where 7 = (j, a,u), and 
a — {b,a) gives the transition | a) — s-j b) that destroys one electron, i.e. with | 6) having one 
electron less than | a) as in Ref. 15l and it is often convenient to abbreviate (7, r) = £. The index 
u is then introduced to describe the inverse transition , so that A (7) = Xj^a when u — — and 
A(7) — xj ^ when u ^ +, and for the PAM with U = 00 there are only two possible transitions 
a = (0,cr). ^ 

The Grand Canonical Ensemble (GCE) of electrons is employed in this problem, and it is useful 
to introduce 

^ = ^-^ E^L^k.^+^i'aA.-aJ , (3.4) 
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where Va is the number of electrons in state | a) and fi is the chemical potential. As usual Ti, is 
split into 

n = Ho+Hh , (3.5) 

where Hh will be the perturbation Hamiltonian, and the exact and unperturbed averages of any 
operator A will be denoted respectively by < A >T-i and < A >. It is also convenient to denote 
the exact or Heisenberg "r evolution" of any operator A with 

A{t) = exp[7^T] A exp[-HT] (3.6) 

and employ A (r) = exp[7io''"] A exp[— TYot] for the unperturbed case. It is convenient to use the 
expressions 

^j,a = Ej^a ~ Vai^ and ek,o- = E]^,a - A* , (3-7) 

because the energies appear almost always in those combinations. 
The c-electron GF employed in the cumulant expansion are: 

G°,(k,c.) = -^ , (3.8) 

where the subindex c was added to differentiate them from the unperturbed GF of the f-electron. 

GlAu;)^--^ , (3.9) 

where -D^J = (X,,) + (Xoo). 



1. The two Hybridization models 
The general hybridization coefficients are given by (cf. Eq. (2.3) in Ref. ||) 

y,-k,. = -i=14(k)exp(ik.Rj) . (3.10) 
The hybridization is purely l ocal when Vo-fk) = V° is independent of k, because using the 



Wannier transformation (cf. Eq.(p.2|)) in Eq. (3^) it follows that 



The Hh that corresponds to n.n. hybridization is 

where the vectors 5 give the position of the d n.n. sites of the origin, and the corresponding VJj(k) 
is immediately obtained: 

V 

K(k) = 2 V, ^cos(k„ a) = — ^ Ek,. . (3.13) 

a=l 
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A. The chain approximation 



In figure |^ are shown some of the infinite diagrams that contribute to the exact GF 
(^Xj^aiT) Xj> a'^ \ . Tlic full circlcs (f- vertices) correspond to the cumulants of the f-clcctrons 



and the empty ones (c-vertices ) to those of the c-electrons. Each line reaching a vertex is associated 
to one of the X operators of the cumulant, and the free lines (i.e. those that do not join an empty 
circle) correspond to the external X operators anaearing in the exact GF. An explicit definition 
of the cumulants can be fopHd_in the referencesD'Sil3,and they can be calculated by employing a 
generalized Wick's theorenJljtZl. 

The first diagram in figure |^a corresponds to the simplest free propagator (^{Xj^aiT) Xji,a')^, 

and the second diagram in that figure has an empty circle that corresponds to the conduction 

electron cumulant, equal to the free propagator ^ ^Cfecr(r) ^. The interaction is represented 

by the lines (edges) joining two vertices and, because of the structure of the hybridization, they 
always join a c-vertex to an f-vertex; the number of edges in a diagram gives its order in the 
perturbation expansion. 

Cumulants containing statistically independent operators are zero, and those appearing in the 
present formalism (with the hybridization as perturbation) vanish unless they contain only X 
operators at the same site or only C or operators with the same k and a. The only non-zero 
c-cumulants are of second order, because the uncorrelated c-operators satisfy Wick's theorem. On 
the other hand, the f-vertices can have many legs, all corresponding to X operators at the same 
site, like the fourth and sixth order cumulants appearing in the rather more complicated diagram 
shown in figure 

All the infinite diagrams that contribute to the GF of the f-electron with cumulants of at most 
second order are shown in figure 0a, and they define the "chain approximation" (CHA) when all 
the other diagrams are neglected. The corresponding approximation for the GF of the c-electrons 
corresponds to the diagrams of figure |l|b. The diagrams of the CHA usually appear as part of more 
complicated diagrams, and it is then useful to analyze their behaviaiir when (i — > oo. In the CHA, 
the GF is given in frequency and k space by (cf. Eq. (3.10) in Rcf.Ej) 

and it would be useful to transform it back to real space to s how the dependence with the distance 
Ri J- between the two sites i and j, as it was done in Eq. ( |2.8| ) for the conduction electrons. Because 
of the lattice translational invariance, it is enough to use the distance Rj of the site j to the origin, 
and characterize this site with an n = (ni, . . . ,nd), (cf. Sec. ||). The set of indices r with 7^ 
will be denoted with {r}j, while s{j)=^fns is the minimum number of n.n. jumps necessary to 
go from the origin to the site j. In both the local and the n.n. hybridization models, the G/_o-(k, uj) 
depends on k only through the Eko- = i'ko- — fJ-, and for any given site j one can write 

Ek, = E,{{r},,k)^0 cos(aks) , (3.15) 

se{r}j 



where 



E,{{r}j,k) ^ -0 cos(aks) (3.16) 



and 6 = t/^/2d. Substituting this relation in Eq. ( [3. 1-1 ) one can expand it in a power series 



of all the cos(a fcg) with s e {r}j and then transform back to real space to obtain Gf^a(Rj,^j- 
Employing the relation 

cos™(fcg a) exp(insksa) — for m < Ug , (3-17) 
which is a consequence of the cyclic boundary conditions, it follows that 
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G/,<,(Rj,w) = O (2d)-'*(^')/2 . (3.18) 

The GF for the conduction electrons in the CHA, corresponding to the diagrams in Fig. |l|b, is 
given by: 

Gc,.(k,c.) = - "^'"""/"Iv^u^ 12 no ' (3-19) 

[luj - ea){iuj - £k<T)- I V(k) p 

and employing the same derivation used above for the f-electrons it follows that GcaCRj,^) = 
0|(2d)-^(^')/2|. The same relation is obtained in the absence of hybridization, as stated at the end 
of Section Q for G° (R„,m, w,.)- 

To close this section let us emphasize that the present expansion employs as perturbation, 
and that the exact solution of the conduction problem in the absence of hybridization is included 
in the zeroth order Hamiltonian. The contribution to the GF joining two sites separated by s{j) 
n.n. jumps, that is of order | (2(i)~'*'^^^/^ |, includes then contdbutions of any order in the He- It 
is because of this difference that the derivation of MetznerO for the Hubbard model had to be 
modified for the present problem. 

All the contributions of He would disappear in the case of a band with zero width, but the 
electronic wave functions would still be extended for the model with n.n. hybridization. 



IV. THE CONTRIBUTION OF DIAGRAMS FOR INFINITE DIMENSION 

In the present section it will be shown that the only diagrams that remain in the limit d —f oo 
are those that are topologically "fully two particle reducible" (f.t.p.r.). These diagrams-are defined 
as those in which any pair of vertices can he-, separated by cutting one or two edgesEa, and are 
"topologies constructed by linking polygons"t3. Two points should be made here: the first is that 
two different polygons (also called loops or rings in what follows) can have at most one vertex in 
common in the f.t.px. topology (see Fig ^,b). The second point, already stressed by Metzner for 
the Hubbard modeljlj is that the topology of a diagram can be different from those of its possible 
embeddings in the Lattice, because in the cumulant expansions there is no excluded site restriction 
in the lattice sumsEl, and two different vertices of a diagram can occupy the same site (see Fig. 
|2|3,d). The property stated above refers to the topology of the embeddings, and in the sum of 
contributions of diagrams that are not f.t.p.r. there may be terms that give non zero contributions 
because they correspond to a f.t.p.r. topology of the embedding: the diagram in Fig. 2c could 
contribute when j = i, because the topology of its embedding, given by Fig. 2d, is f.t.p.r. 

To prove this property it is necessary to modify the derivation given by Metzner for the Hubbard 



model, for the same reasons given in Section III A for the GF dependence on s(i). The proof of this 



property is given in the following two subsections, and it can be summarized as follows. Consider 
first the diagram in Fig. ||a, that corresponds to the sum of loops with all possible lengths. All the 
embeddings of these diagrams are f.t.p.r., and they give a non-zero contribution to the calculation 
of the free energy. As a second step, one can consider diagrams like those in Figs. ^,d, obtained 
by linking loops of any length at different f-vetjtices. One could obtain all these diagrams by a 



procedure similar to the vertex renormalizatio with the difference that only insertions of a 

rather special type are considered. All these diagrams, as well as their embeddings, are f.t.p.r., and 
a finite contribution of all of them is expected. The next step is to consider diagrams that are not 
f.t.p.r., like the one shown in Fig that can be obtained by joining two loops (like that in Fig. 
|3^) at two different f-vertices. If the family in Fig. ^ gives a finite contribution, one can show that 
all the contributions to the diagram in Fig. ^ will vanish when the topology of the embedding 
coincides with that of the diagram, i.e. when the embedding is not f.t.p.r., as it would happen 
when i ^ j. In the special case of i = j, (cf. Figs. || c,d for a special case of this situation) the 
embedding is f.t.p.r., and the diagram contribution does not necessarily vanish. Of the two sums 
over the lattice sites i and j, only one over i — j remains in the limit d ^ oo, and this simplifies 
the calculation of this diagram. In particular, fhs, restriction that in reciprocal space there must 
be conservation of k at each vertex is removeclQllJ, leaving only a single conservation of k for all 
the edges joining the collapsed vertices: this makes the calculation of the diagram much simpler 
in the limit c? — > oo. 
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A. The f.t.p.r. diagrams 



To prove that the contribution of a f.t.p.r. diagram can be finite, let us consider the family of 
diagrams represented in Fig. ^a. Their contribution does not vanish and can be expressed as 

J Ml j M2 ^°(ji;^i,^2) M°(ji;^i,£2) (4.1) 

h 

where 

M°{n:hA2) = ((X (£1)^(^2))+), (4.2) 

is a local cumulan1@ and / Ms = X^a jo ■ The abbreviations s = £s and J ds = J dig will 
be used when there is no possibility of confusion. The symbol ^'^(ji; ^1, ^2) = <S'2(ji;l,2) corre- 
sponds closely to the "self-fieWB S*™ ( ji ; ^1 , ^2 , ■ • ■ , tm) = Sm{ji', 1,2,..., m), that was employed 
in Ref. ^ to renormalize vertices, but it gives only "insertions" obtained from simple loops of any 
length, as represented by the diagrams in Fig by it definition, all the Is correspond to the 
same site, indicated by ji. The notation M2{ji;ii,£2) = -^^2(^1! l^-^) has the same meaning, and 
reflects the fact that this cumulant is zero unless ji = j2 because X{£i) and -^(^2) are statistically 
independent in the unperturbed system when ji ^ ji- Note that S'2 (ji; ^i, £2) depends explicitly 
on the parameters (where ii represents ji,Ui,ai — {0,ai) for z = 1,2) through the two 

hybridization constants Vj.k.cr associated to the edges joining the insertion vertex. 

To consider the addition of a simple loop to any f- vertex of a diagram, it is convenient to consider 
first the simplest possible case, shown in Fig. its contribution is obtained by substituting the 



cumulant M^{ji] 1,2) in Eq. 4.1 by 



JdsJdA M°(j-i; 1, 2, 3, 4) 52°(ji; 3, 4) . (4.3) 

When the loop is added to a site ji that is already joined by n edges, the corresponding cumulant 
M^(ji; 1, 2, . . . , m) suffers a similar substitution. Repeated application of this procedure at all f 
vertices gives all the possible f.t.p.r diagrams, and as the cumulants are independent of the lattice 
dimension d and of the site j, this procedure should not affect the order of the contribution with 
respect to d. The simple loops of Fig. ^ give a finite contribution, so that the contribution of any 

f.t.p.r diagrams is then of O (^\df^. 

A similar type of procedure can be applied to the CHA diagrams of the one-particle GF, by suc- 
cessive decoration of the f- vertices with any number of insertions corresponding to the >5'2 (ji; ^1, ^2) 
discussed above. It is then clear that these diagrams will still be f.t.p.r., and that the corresponding 
GF joining two sites separated by s n.n. jumps will be of order O |(2(i)~*/^| , as it was shown in 
Section Iml for the CHA. 



B. The diagrams that are not f.t.p.r. 

To analyze the diagrams that are not f.t.p.r., consider any one of them as a "mother" diagram, 
and split it into two or more f.t.p.r linked "daughter" diagrams without any edges in common but 
such that any of them has at least two vertices in common with another daughter diagram, as well 
as two edges arriving at each of the common vertices. It is clear that any pair of vertices that are 
common to two daughter diagrams can not be separated in the mother diagram by cutting one or 
two edges, and one says that they are not "two particle reducible" (t.p.r.). 

From all the daughter diagrams choose one as a vacuum diagram and transform the remaining 
ones into GF by adding external lines of the Bose type (as indicated in Appendix ^ to all the 
vertices that each of them had in common with any other daughter diagram in the mother diagram. 
Three examples of this procedure are given in Fig ^. 

Consider now the real space calculation of the contribution of the daughter diagrams. In Ap- 
pendix ^ it is shown that by joining one daughter GF to another diagram and fixing the position 
of the common vertices, one obtains the order of the contribution of the resulting diagram as the 
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product of those of the building diagrams. In Appendix ^ it is shown that unless all the external 
f-vertices of one daughter GF coincide at the same lattice site, the corresponding contribution is 
of O 1(0)^1 with p > 1, where 9 = t/^/2d. Taking now the only daughter vacuum diagram and 
adding successively all the other daughter GF diagrams, the total contribution for given positions 
of all the GF external vertices will then be O \ {9Y\ with p > 1, unless all the external vertices 
of each daughter GF coincide at the same lattice site, that can be different for different daughter 
GF. A typical example of this type of diagram is given in Fig |^, while in Fig ^jc it is shown 
the corresponding f.t.p.r. topology of the embedding when the vertices that are not t.p.r. in the 
diagram are at the same lattice site. 

It is then clear that in the sum over the vertices that are common to all the daughter diagrams, 
each of the terms will be at best of O and would vanish in the limit d ^ oo, unless all the 
external vertices of each daughter diagram are at the same lattice site, that can be different for 
different daughter diagrams. One can conclude that the contribution of the mother diagram would 
vanish when any pair of vertices i and j that are not t.p.r. occupy different lattice sites, because in 
that case, it is always possible to separate from the mother diagram a daughter GF that has i and 
j in common with the rest of the diagram. In that case, only a sum over i — j should remain from 
the unrestricted sum over i and j, and those two vertices then become t.p.r. in the embedding of 
the diagram (this is illustrated in Figs. ^DjC). This reduction of the terms that contribute to the 
lattice summations has been already described|,|-and called "collapse of vertices", in the study of 
the U-perturbation theory in high dimcnsionalfl'ta: this property will also be given the same name 
in the present work. 

One important point to notice, is that the collapse of vertices occurs only in the embedding of 
the diagram, and that the rules for calculating the contribution should be applied to the original 
diagram (i.e. it would be the diagram of Fig. |^ in the examples given) and not to the collapsed 
diagram of the embedding (i.e. the diagram in Fig. |^). 

As in the collapse of i and j the two independent lattice summations over the vertices i and 
j are replaced by a single lattice summation over i = j, it is easy to see from the derivation of 
the contribution rules in reciprocal space (cf. Section III D in Ref. ||) that the two independent 
momentum conservation at each of the two collapsed vertices becomes a single conservation of the 
wave vectors corresponding to all the internal edges joining them. When all the pair of vertices 
that are not t.p.r. have been collapsed, one can see that at every vertex of the embedding of 
a vacuum or one particle GF diagram, all the edges can be arranged in pairs with momentum 
kg and -kg respectively (cf. Fig. ||c as a typical example) It then follows that the momentum 
conservation at the collapsed vertices are automatically satisfied, and most of the usual restrains 
on the momentum integration disappear in the limit d ^ oo. 



V. SUMMARY AND CONCLUSIONS 

The cumulant expansion of the PAlVli was considered in the limit of infinite spatial dimensions 
(d oo) for two types of hybridization models: the purely "local" hybridization and the "nearest 
neighbor" (n.n.) one. As the systems usually described by the PAM have a bandwidth much 
larger than the hybridization energies, the unperturbed Hamiltonian is chosen to include the exact 
solution of the conduction band electrons inptte absence of hybridization. It was then necessary to 
modify the derivation employed by MetznerEJ for the cumulant expansion of the Hubbard model, 
in which the hopping Hamiltonian is used as perturbation. The basic result presented here is that 
in spite of this change the PAM shows, for the two hybridization models considered, the same type 
of simplifications that occur in the diagrammatic cumulant expansion of the Hubbard model when 
d — > oo. 

Only the linked diagrams contribute to the general cumulant expansion, and there is a vertex 
collapse when d —^ oo for those diagrams that are not f.t.p.r. A diagram is not f.t.p.r. when at 
least two vertices i and j can not be separated by cutting only one or two edges of the diagram, and 
their collapse means that from the two independent summations over i and j that are necessary to 
calculate their contribution, only a summation over i = j remains. The topology of the embedding 
is different from that of the diagram itself when i = j, and the collapse is repeated until the 
topology of the embedding becomes f.t.p.r. One important point to notice is that, after the collapse 
of vertices, the rules for calculating the contribution should be applied to the original diagram and 
not to the collapsed diagram of the embedding. 
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When the two independent summations over i and j of the diagram become a single summation 
over i = j because the two vertices i and j collapse, the two separate conservations of k at these 
vertices become a single conservation of the k corresponding to all the edges joining both i and 
j. As a consequence, most of the usual restrains on the momentum integration disappear in the 
limit c? — > oo, and the calculation of many diagrams is very much simplified by this change. One 
should note that the vertex collapse does not alter the calculation over frequencies, which keep 
their conservation at all the diagram vertices, even when cL=> oo. 

Employing the cumulant expansion, Metzner has shownE£l that the single-particle properties of 
the Hubbard model in the limit c? — > oo, can be described as that of independent electrons "hopping 
between dressed atoms characterized by an effective Green's function" . A similar derivation can be 
employed for the PAM, and the exact one electron GF is given by the family of diagrams in Fig. |l|a, 
but using an effective cumulant M^^/ (oj) for the f-elcctron vertices instead of the bare one AI^ ai^)- 
The effective cumulant M2^J (to) is given by the contribution of all the diagrams of G'CT(Rn,m — 0, uj) 
that can not be separated by cutting a single edge (usually called "irreducible diagrams"), where 
^^^•(Rn.m = 0, bj) is the exact GF of the f-electrons in the real space representation for n = m. This 
property is only valid in the limit d — > oo, so that all the diagrams that are not f.t.p.r have their 
vertices collapsed until the associated embeddings are f.t.p.r. The exact GF can then be written 

Gf<,(k,t^) =M|{.-^(w) ^- (5.1) 

A practical difficulty in the study of correlated electron systems with cumulant expansions, is that 
the higher order cumulants rapidly become very laborious to calculate. To consider in some way the 
higher order cumulants, we are studying the substitution of M^^J {ui) by an approximate quantity 
M2ct(w)j derived from an exactly soluble model. To this purpose we use the same Anderson 
periodic model but in the atomic limitll3, i.e. when the hopping of the conduction electrons is 
eliminated by taking a conduction band of zeroth width. This attempt will be discussed in another 
publication, and the present work provides an essential guidance to that study, by showing that all 
the cumulant diagrams present in M2^J {uj) are also present in the approximate M2*^{lu). Although 
the hopping is missing from the Mf '^.(wl. the exact solution of the conduction band shall appear 



in the "hopping between dressed atoms"L3 through the G°^(fc,a;) in Eq. (5.1 ) 
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APPENDIX A: SEPARATION OF VERTICES 

In Section |^ it is shown that the contribution of a "mother" vacuum diagram that is not f.t.p.r. 
can be estimated from those of several f.t.p.r. "daughter" diagram that are derived from the 
original one. The daughter diagrams are obtained by separating the mother diagram into several 
f.t.p.r. subdiagrams: a vacuum diagram plus several GF diagrams obtained by adding external 
lines of the Bose type to the vertices that these separated diagrams have in common in the mother 
diagram. To show the procedure, it is convenient to analyze a simple case, like that of Fig. 
and consider only a single vertex in real space, as shown in Fig assuming first that it does not 
belong to the vacuum daughter diagram. Writing only that part P of the total contribution to the 
mother diagram that have parameters connected to the particular vertex under consideration, one 
obtains 

P = ^ y dl y d2 y d3 y d4 ^ ^ exp [-i (ki - ka + kg - k4) .Rj] 

j kik2 k3k4 

xy(l)F*(2)T/(3)F*(4)G0(l)G°(2)G°(3)G^(4)M0(j,l,2,3,4) (Al) 
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where the G°(s) are the c-electron GF, and the other abbreviations are those introduced in Section 



IV , The time or frequency dependence has been left out because it only plays a trivial role in this 



proof. The corresponding contribution Pa to one daughter diagram for fixed values li,£2 of the 
internal lines joining the vertex, and with an added external line of momentum Ka and parameters 
L is 

Pa(Ka, (a, ^1. ^2) = ^ + '^j] 



kik2 J 



xMl{j,a,l,2)Vil)V*i2)G°il)G°{2) , (A2) 

and the contribution Pfc(Kb, ^h, ^3, £4) of the other daughter diagram is trivially obtained by re- 
placing a, 1, 2 by b, 3, 4. 

Transforming Pa(Ka, -^a, ^1, ^2) to real space, the contribution associated to a site Rj is given by 



Pa{Rj,£aJl,£2) = ^^exp {iKa.Rj) PaiKaJa) 



= J2 exp [-Z (ki - k2) .R,] 1-(1)T/*(2)G0(1)G°(2)M3"(j, a, 1, 2) (A3) 

kik2 

and a similar expression Pf,(Rj, £b) is obtained for the other daughter diagram. It is clear that the 
only difference between the product Pa{'R.j,£a,£i,£2) -Pfc(Rj , ^6, ^3, ^4) and the term corresponding 
to Rj and fixed ^1,^2,^3,^4 in Eq. Al is given by the cumulants M", which are all independent of 
both Rj and of the space dimension d. It then follows that 

P^ ldlld2ld3lM ^^^°(^'1'2'3,4) 



X 

3 



Ml{j,a,l,2)MlU,b,3A) 
^Pa(R,,4,^i,^2) P6(Rj,4,4,4) (A4) 



gives the dependence of P with d. 

When one of the two daughter diagrams is the chosen as the vacuum one, it is necessary to 
consider the corresponding contribution 

P{Rj,£u£2) = ^ expH(ki-k2).R,]y(l)^^*(2)G:?(l)G°(2)M0(j,l,2) , (A5) 

kik2 

that was obtained replacing M^{j, a, 1,2) by M^j, 1,2) in Eq. ([A^): this is the quantity that 
appe ars in the corresponding vertex of the vacuum daughter diagram. In this case, instead of Eq. 
( A4) one has the relation 

P^ Idlld2ld3ld, ^^°(^'1'2'3,4) 



X 

3 



M0(j,l,2)M0(j,6,3,4) 
^P(R,,£i,^2) n(Rj,4,^3,4) . (A6) 



This are the relations employed in Section IV to prove that only embeddings that are f.t.p.r. give 
a non zero contribution, thus leading to the collapse of vertices when they can not be separated 
by cutting at most two edges. 



APPENDIX B: DIAGRAMS THAT ARE NOT F.T.P.R. 



In the procedure discussed in Section IV a "mother" vacuum diagram that is not f.t.p.r. is 
separated into several subdiagrams that have that property. As discussed before, one of the 
daughter diagrams is chosen as a vacuum diagram, and the remaining ones are GF diagrams 
obtained by adding external lines of the Bose type to the vertices that these GF diagrams share 
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with other daughter diagrams in the mother diagram. Let us consider one of the GF daughter 
diagrams in the reciprocal space, with momenta K„, indices and transitions a^, {v = 1,2, . . . ,n) 
assigned to the n> 2 external vertices, where the a-u correspond to pair of states with equal number 
of electrons, i.e. to operators Xa of the Bose type. All these external parameters will be indicated 
by {Kv,£v}, while {i}int will be used for the set of £ associated to the internal edges joining all 
the external vertices. The contribution corresponding to fixed {£}int of a split diagram is denoted 
by F ({K„, £1,} , {£}i„t) /iV"/^, where the factor 1/^/N, associated to each external line joining an 
f-vertex (cf. Rule 3.7 f in Ref. ^), has been explicitly written. To transform from K^; to position 
variables one should calculate 



F{{Il,,},{£}) 



[f{{K,,£,},{£}^„^)/N"/' 

{K„} 



exp 



v=l 



(Bl) 



where {£} denotes all the £, both internal and external. From the translational invariance of the 
system (or from wave vector conservation at all the vertices) it follows that J^v = 0, so that 



F ({R4 , {£}) = I n (^^ E (R. - Rn)]^ I 



(B2) 



In the two hybridization models considered in this work, the external wave vectors appear in 
F {{K^,£^} , {-^Ijnj) only through the delta that gives the wave vector conservation at each external 
vertex. The internal wave vectors k appear through the E'lco- in ^{uj) (i.e. the c-electron GF 
of Eq.dU)) and also in the K(k) = - (V^/t) Ek,<T (cf. Eq. (|3.ic[ )') for the n.n. hybridization 
model. When the wave vector conservation at all vertices is considered explicitly, the number of 
summations over internal k is reduced, but the arguments of the iJk.a- become linear combinations 
of both the remaining external and internal k. To illustrate this result, consider the diagram in 
Fig. 0: applying momentum conservation one can write 



({K„ 4, Ka, 4}, {^1,^2,4,4}) = 
J2 Ml{b, 1, 2) Ml{a, 3, 4) \V4k^)f |F.(ki - Ka)!' G° (ki) G° (ki - Ka), 



(B3) 



ki 



where the frequencies are not explicitly written because they don't play any role in this analysis. 
There is only a single summation over ki , because from the sum over k2 only k2 = ki — Ka 
remains. 



To study the general case, it is convenient to write (cf. Eq. (2.4)) 



Ei^^a = -Bo-(k) = - cos (a kg) 



s=l 



(B4) 



so that F ({K„, 4} , {^}i„t) can be expanded in series of the 9s, putting 9^ = 9 = t/V2d at 
the end of the derivation. To analyze the dependence of F {{Ri,} , {£}) with 9, it is conve- 
nient to concentrate first on a given external vertex v, and write = Q = (Qi, • • ■ iQd) and 
— R„ = R = {Ri,...,Rd)- The Eq. (B2) shows that for the given external vertex v, 
there is a X^g^ ^xp (z QsRs) apphed to F {{K.v,£i,} ,{£}^^^) for each dimension s = of 
the space, and when for a given s it is i?<; 7^ 0, all the terms independent of 9s in the expansion of 
F{{K^,£^} ,{£}^^J will cancel out (cf. Eq. (p7|)). It then follows that F{{R^},{£}) =0\9p\, 
where p > and p^ is the number of dimensions for which there arc non zero components s for 
at least one of the n — 1 vectors R^, — R„. One can then conclude that unless all the coincide, 
the F ({Rt,} , {^})is at least of O \9\ and it vanishes when d 00. This property was used to prove 
the collapse of vertices discussed in Section [V for diagrams that are not f.t.p.r. 
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FIG. 1. Typical cumulant diagrams for one-particle GF (a) The diagrams of the chain approximation 
(CHA) for the f-electron, represented by the filled square to the right, (b) As (a) but for the c-electrons 
, represented by an empty square, (c) A more complicated diagram, with cumulants of fourth and sixth 
order. 



FIG. 2. Examples of relevant topologies in the limit d ^ oo. a)FuIly two particle reducible (f.t.p.r.): 

any pair of vertices can be separated by cuting two edges. Note that all pairs of loops have at most one 
vertex in common, b) The diagram is not f.t.p.r. c) The topology of the diagram and that of the embedding 
are the same when i =^ j. d) The embedding topology that corresponds to the diagram in c) when i = j. 



FIG. 3. a) The sum of all the simple loops of any length, b) Diagrams of an insertion with contribution 
S^iji'-i^iT^i), obtained by fixing the site ji and the parameters £i,£2 of one conduction vertex (that is 
considered external) of the diagrams in a). The full circle in these vertices, corresponding to the cumulant 
M2 (ji ; £i, ^2), is removed from the diagrams and is replaced by 1 in S2{ji', £1, £2). c) A f.t.p.r diagram 
obtained by joining two simple loops at an f-vertex. d) Family of four linked loops that give a ft.p.r 
diagram, e) Family of diagrams obtained by linking two loops at two different f-vertices: it is not ft.p.r. 



FIG. 4. Three examples of separation of a mother vacuum diagram into several daughter diagrams, a) 
Simplest case of separation with two common vertices. b),c) Separation in three daughter diagrams with 
a total of four shared vertices. 



FIG. 5. a) A family of f.t.p.r diagrams that give a finite contribution for large d. b)Family of diagrams 
obtained by joining pairs of f.t.p.r diagrams in two common vertices; the diagrams are not f.t.p.r. c) The 
f.t.p.r topology of the embedding of b) when i = j. Note that the diagram contribution is calculated with 
the diagram b) and not with c), but replacing the two independent summations over i and j by a single 
one over i = j. 



FIG. 6. Separation of a vertex with four edges into a pair of vertices with two edges each. 
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FIG. 7. A simple f.t.p.r diagram. The conservation of moment k at the vertex b gives ki — k2 — Kb = 
for Ub = —1 and that at vertex a gives k2 — ki + Ka = for = +1, so that Ka = Kb as should be 
expected from the translational invariance of the system. Only one summation over ki remains because 

k2 = ki - Ka. 
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